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Abstract. 

Transverse elastic waves behave differently in nonlinear isotropic and anisotropic media. While 
^ ■ in the former the quadratically nonlinear coupling in the evolution equations for wave amplitudes is 

not possible, such a coupling may occur for certain directions in anisotropic materials. We identify 
the expression responsible for the coupling and we derive coupled canonical evolution equations for 
transverse wave amplitudes in the case of the two-fold and three-fold symmetry acoustic axes. We 
illustrate our considerations by examples for a cubic crystal. 
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C^ 1. INTRODUCTION 

It is known that transverse elastic waves cannot interact with each other on the quadrati- 
cally nonlinear level in an isotropic material [5]. On the other hand it is also known that 
>■ ■ in anisotropic materials this is not the case [1]. There are certain directions in crystals 

Jf[^ . for which quadratically nonlinear interactions do occur between (quasi-) shear elastic 

\Q \ waves. This manifests itself by the presence of couplings in the evolution equations for 

■^ I pairs of wave amplitudes. We show that the two-fold and three-fold symmetry acous- 

t:j- ■ tic axes are examples of these special directions along which such a coupling takes 

O \ place in crystals. We identify the quantities which when nonzero are responsible for 

Q \ the occurrence of the coupling between nonlinear elastic shear waves. Then we use the 

perturbation method of weakly nonlinear geometric acoustics and we derive coupled 
canonical evolution equations for shear wave amplitudes of plane waves propagating 
rS \ along two-fold and three-fold symmetry acoustic axes. We illustrate our considerations 

c^ ■ by examples for cubic crystals. 

2. PLANE WAVE EQUATIONS 

Let us consider the system of equations describing plane waves propagating in n direc- 
tion in an arbitrary nonlinear anisotropic hyperelastic material in which both geometrical 
and physical nonlinearities are taken into account. This system of equations can be writ- 
ten as follows [1, 2, 3, 4]: 

w,f + A(w,n)w,x = 0, (1) 
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with 



w 



v{x,t) 
m{x,t) 



*/ \ I -f B(m,ii) , ,-, 

A(w,n) = -( J P\ ], (2) 



where 

B = A + ^m + --. (3) 

Here Po is a constant density \ v is a velocity, m is a 1-D strain, n is the direction of the 
plane wave propagation, I is a 3 x 3 identity matrix and B is a 3 x 3 symmetric matrix, 
which we restrict to consist of zero and first order terms in strains A and ^, respectively. 
Alternatively, in components, 

Bac = Cabcdnbnd + CahcdefUbUcinf nig -\ . (4) 

where Cabcd and Cabcdef ^e second and third order elastic constants and Ua and nia are 
components of the vectors n and m, respectively. By the standard reduction of indices 
we represent the second order constants with two indices and third order constants with 
three indices. 

We assume that the Christoffel tensor has the form 

A= aiki(g)ki + a2k2(S)k2 + a3k3(g)k3, (5) 

where Uj > and {ki,k2,k3} is an orthonormal triad of vectors. The six eigenvalues of 
A(0) are the characteristic speeds of a pair of longitudinal and two pairs of shear waves. 



(6) 
(7) 
(8) 



r2/-i 



•2,-1 



Ai 


= -V^i = -A2, 


A3 


= — v^0f2 = — A4, 


A5 


= -^"3 = -h- 


left eigenvectors are, respectively. 


/^-A2,_ikA 




1 (-^2jlM 

2 V k; ; ' 


1,^. = 1M2> 
' 2^ k, 



7 = 1,2,3, (9) 



We have 1, ■ r^ = 5,j. We are interested in cases when the matrix A(0) has coincident 
eigenvalues corresponding to the speeds of shear waves. Such waves will be called 
degenerate. Hence suppose that a2 = OC3, so A.v = Xs+2 for 5 = 3,4. In such a case n 
is an acoustic axis (see [7], or [10]). 



which, in what follows, we assume to be equal to 1 for simplicity. 
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3. COUPLED EVOLUTION EQUATIONS 

Let us consider the initial- value problem for the system of plane waves elastodynamics 
(1), (2) with perturbed initial data which are of compact support (£ is a small parameter): 

w^j+A(w^,n)w^^ = 0, 

(10) 

v/''\t=Q = ey/i{x,x/£). 

We apply the asymptotic method of weakly nonlinear geometric acoustics. The method 
combines the use of a multiscale perturbation with a small amplitude, high frequency 
asymptotics [1]. Assume that the asymptotic expansion is of the form 



y/{t,x)=£\Os{t,x,r\)rs + Os+2{t,x,r\)rs+2\ +0(e^), 5 = 3,4, (11) 

with 9 = e^^{x — A,t), A = A.v = As+2, and where Os, (7^+2 are the unknown amplitudes. 
Inserting (1 1) into (10) and applying the method of multiple-scale asymptotics (see [4]), 
we obtain the following coupled evolution equations for the amplitudes a^ and Os+2'- 

(c7,),, + A,(c7,),x + r,c7,(c7,),e +rf ^(cT.s-cJ,+2),e +n+2C^.v+2(cy.s-+2),e = 0, 

(o".s-+2).f + A,v(asO,x + r;::+2o-,(c7,),e +n+2(o".s-«^.v+2),e +r,+2ty,v+2(o-,+2),e = 0. 

(12) 
Here we denote F;^ ^ = F^, and we use the relations which hold for elastodynamics 

T^s+2 — -r^s r^.s+2 -pi- — r^.s+2 -pi- n ^\ 

^s =^s,s+2 — ^s,s 7 -^ .s+2 = -^ .s,.s+2 ~ -^ .v+2,i-+2' V^-'-' 

where the interaction coefficients are defined 

r;,, = lr(VwA(w)rpr,)|^^0- (14) 

3.1. Two-fold symmetry acoustic axis 

In the case of a two-fold symmetry acoustic axis, the evolution equations further 
simplify and it turns out (see [4]) that only two coefficients Fy and T'Ij_2 characterize 
the nonlinear terms in the coupled evolution equations (12) for a pair of^ amplitudes of 
(quasi-)shear waves, in that case. 



3.2. Three-fold symmetry acoustic axis 

In the case of a three-fold symmetry acoustic axis only one coefficient Fs charac- 
terizes the nonlinear terms in the evolution equations for a pair of amplitudes of plane 
degenerate shear waves, and the coupled evolution equations are as follows: 

((y,v)j + Av(cy,v),x + F, [(J,(C7,),0 - (y,+2{Os+2),e] = 0, 
(0',-+2),f + As-(0-,-+2),x - F,-(0-,- as-+2),0 = 0. 
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4. EXAMPLES 

In the examples below we consider a cubic crystal of class m3m in which the strain 
energy W is characterized by three second order and six third order elastic constants 
(see [1, 2, 3, 4]) 

W = W{cn,Cl2,C44,Cn\,Cn2,C\44,C\23,Cl66,C456)- 

Example 1. Here we present the case where the propagation direction of the plane 
wave is a two-fold symmetry axis n = 4= [1 1 0] which is not an acoustic axis. The speeds 
of shear waves are (see [1, 2, 3, 4]) 



C11-C12 



As 



-A, 



hence A3 7^ A4 7^ A5 7^ Ag, and Ts = for 5 = 3, 4, 5, 6. It turns out that there is no coupling 
and moreover one can show that decoupled evolution equations with cubic nonlinearity 
describe the propagation of shear waves in this case (see [1]). 

Example 2. Consider now the case of a three-fold symmetry axis n = 4= [1 1 1] which 

is the acoustic axis. In this case we obtain (see [1, 2, 3, 4]) coupled equations (15) with 



A3 = A5 



Cll -C12 + C44 



— Aa = — Afi. 



r3 = 



18V2 



Cm +2C123-2C456-3(C112-Ci44+C 



166 J 



A3 



= -r4. 
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